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1.0  INTRODUCTION 


Applied  Field 

Fixed  Magnets 

Fixed  Magnets 

Cooling 

- 

Liquid 

Flow 

- 

1/4  gpm 

Pressure  Drop 

- 

10  psi 
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Initially,  a  four  port  differential  phase  shift  circulator 
using  nonreciprocal  phase  shifter  in  a  rectangular  waveguide 
structure  was  studied.  Nonreciprocal  phase  shifter  design 
with  H-plane  ferrite  and  E-plane  dielectric  in  a  rectangular 
waveguide  displayed  limitations  in  terms  of  broadband  performance 
With  optimum  phase  shifter  geometry,  90%  of  required  bandwidth 
could  be  achieved.  ^However  matching  section  of  -phase  shifter  and 
transition  from  rectangular  to  ridge  waveguide  would  further 
limit  the  bandwidth.  Therefore,  the  phase  shifter  design 
using  rectangular  waveguide  was  abandoned  and  nonreciprocal 
phase  shifters  in  double  ridge  waveguide  structures  were  used  in 
the  circulator  design. 

-  During  the  course  of  the  circulator  development  program 
analytical  studies  on  dielectric  and  ferrite  loaded  double 
ridge  waveguide  were  performed.  Design  effort  was  placed  on 
the  frequency  range  between  8  to  16  GHz.  Raytheon  circulator 
design  uses  two  magic  tees  and  two  non- identical  nonreciprocal 
phase  shifters  in  standard  V'JRD750D24  double  ridge  waveguide. 


2.0 


CIRCULATOR  CONFIGURATIONS 


Differential  phase  shift  circulators  can  be  built  in 
several  configurations.  Four  configurations  are  shown  in  Figure  1 
(a),  (b),  (c),  and  (d).  All  four  have  the  sane  principle  of 
operations.  Power  entering  Port  1  emerges  from  Port  2  only,  and 
so  on  until  power  entering  Port  4  emerges  from  Port  1  only. 

The  device  in  Figure  1(a)  is  the  most  common  configuration. 
This  option  is  probably  the  most  favorable  if  a  3  dB  quadrature 
hybrid  operating  over  the  required  frequency  range  with  acceptable 
unbalance  is  available.  Unfortunately  a  3  dB  quadrature  hybrid 
at  the  desired  frequency  range  with  suitable  performance  is 
not  currently  available.  The  design  indicated  in  Figure  1(b)  uses 
input  and  output  magic  tees,  separated  by  a  0°/180°  nonreciprocal 
phase  shifter  in  parallel  with  a  dielectrically  loaded  reciprocal 
tracking  phase  shifter.  The  design  of  0°/180O  nonreciprocal  phase 
shifters  requires  longer  ferrite  slab,  hence,  causing  higher 
insertion  loss  of  the  phase  shifter.  The  circulator  design 
using  Figure  1(b)  causes  an  amplitude  unbalance  because  the 
insertion  loss  of  the  ferrite  is  higher  than  that  of  the  dielectric 
The  design  in  Figure  1(c)  uses  two  magic  tees,  two  identical  0°/90° 
nonreciprocal  phase  shifters,  and  two  dielectrically  loaded 
reciprocal  phase  shifters  which  provide  a  fixed  reciprocal 
phase  difference  of  90°  between  the  two  phase  shifters 
over  the  operating  frequency  range.  For  wideband  design,  the 
reciprocal  phase  shifters  imposes  serious  limitations  on  the 
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bandwidth  because  a  phase  shift  of  approximately  90°  cannot  be 
maintained  across  a  wide  range  of  frequencies.  The  approach 
using  Figure  1(d)  uses  two  magic  tees  and  two  non-identical 
nonreciprocal  phase  shifters.  Phase  tracking  requirement 
over  wide  bandwidth  can  be  net  more  easily  using  this  technique 
This  approach  v/as  used  in  our  design.  Scattering  matrix 
analysis  for  this  configuration  is  shown  in  the  next  section. 


2 . 1  Scattering  Matrix  Analysis  of  Circulator 

Dynamic  behavior  of  a  four  port  differential  phase 
shift  circulator  can  be  described  using  scattering  matrix 
analysis.  Employing  scattering  matrix,  wave  amplitude  emerging 
from  all  four  ports  can  be  described  as  a  function  of  input 
wave  amplitude  and  scattering  parameters  of  all  the  inter¬ 
connected  devices.  Figure  2  shows  four  port  differential 
phase  shift  circulator  configuration  used  for  anlaysis. 

Instead  of  using  a  conventional  S-parameter  notation  such  as  Sj^ 
a  slightly  different  notation  is  used  in  this  analysis.  They 
are  defined  as  follows: 

Xi  -  indicates  wave  amplitude  entering  into  Port  i  of  the 
magic  tee  1. 

Xi  -  indicates  wave  amplitude  leaving  from  port  i  of  the 
magic  tee  1. 

a^j  -  indicates  scattering  parameter  of  the  magic  tee  1. 

Yj^  -  indicates  wave  amplitude  entering  into  Port  i  of 
the  magic  tee  2. 

Yi  -  indicates  wave  amplitude  leaving  from  Port  i  of 
the  magic  tee  2. 

bj^j  -  indicates  scattering  parameter  of  the  magic  tee  2. 


Ui  - 


Ui  - 


Cij 


Vi  - 


Vi  - 


'i: 


indicates  wave  amplitude  entering  into  Port  i  of 
the  phase  shifter  1. 

indicates  wave  amplitude  leaving  from  Port  i  of 
the  phase  shifter  1. 

indicates  scattering  parameter  of  the  phase  shifter 

1. 

indicates  wave  amplitude  entering  into  Port  i  of 
the  phase  shifter  2. 

indicates  wave  amplitude  leaving  from  port  i  of 
the  phase  shifter  2. 

indicates  scattering  parameter  of  the  phase  shifter 

2. 


All  of  the  parameters  defined  above  are  complex 
quantities  involving  magnitude  and  phase.  According  to  the 
notations  given  above,  each  device  can  be  expressed  as: 


‘ 

r 

^2 

X3 

.  ’'4. 

■  yr 

y2 

y3 

.^4. 

- 

■31I 

1- 

^2 

^2 


■311 

^21 

^31 

341 

'^11 

b21 

^31 

b4i 

^11 

^21 

dll 

^21 


^12 

^22 

332 

342 

t>12 

b22 

t»32 

b42 

C12' 

C22 

<^12 

^22 


313 

323 

333 

343 

^13 

>323 

b33 

b43 
Ui 


^2 

^3 

- 

.^4. 

'Yl' 

Y2 

Y3 

.^4. 

■vi' 

.V2 

314 

324 

334 

344 

bi4 

fc>24 

b34 

b44 


For  phase  shifter  1 


For  phase  shifter  2 


For  magic  tee  1 


For  magic  tee  2 


(1) 


(2; 


(3) 


4) 


And  we  have  the  following  relations: 


X2  =  ui 
X3  =  VI 

(5) 

Y2  =  U2 
Y3  =  V2 

Ui  =  X2 

U  2  =72 

(6) 

Vi  =  X3 

V2  =  73 

VJe  can  identif7  X]^  ,  X4 ,  Y]^,  and  Y4  as  input  wave  amplitude  of 
the  circulator  and  xj^,  X4 ,  yi,  and  74  as  output  wave  amplitude 
of  the  circulator. 

Using  equations  (1)  through  (6),  an7  output  wave 
amplitude  can  be  expressed  as  a  function  of  input  wave  amplitude 
and  scattering  parameters  of  all  the  interconnected  devices. 

For  example,  if  we  consider  Port  1  as  an  input 

port  then: 

X4  =  Yi  =  Y4  =  0  (7) 

and  output  wave  amplitude  can  be  written  as: 

XI  =  aiiXi  +  ai2  (^11X2  +  =1272)  +  ^13  (^11x3  +  di273) 

X4  =  341X3  +  342  (C11X2  +  Cl2y2)  ^43  (^11^3  <^12^3)  (®) 

71  =  bi2  (C21X2  +  C22y2)  +  ^>13  ('^21^3  '^22^3) 

74  =  b42  (C21X2  C22y2)  *=43  (^^21^3  ^22y3) 
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where  X2t  X',  y2»  and  73  can  be  obtained  by  solving  the  following 
simultaneous  equations: 

(l-a22Cii)  -323^11  '322^12  -^23^12  ^2  ^2l'^l 

-332^11  (1-333^11)  '^32^12  "333^12  ^3  ^31^1 

-b22C21  -^>23=^21  (1-022^22)  “^>23^22  ¥2  0 

-t>32C2l  -b33‘321  "'^32=22  ( 1“*^ 33^  22 )  J  |_y  ij  L  - 

Similarly,  we  can  obtain  all  the  output  wave  amplitude  involving 
with  different  input  port. 

If  we  assume  ideal  magic  tee,  scattering  matrix 
of  ideal  magic  tee  can  be  written  as: 


0 

1/v  2 

1//2 

0 

l//'2 

0 

0 

l/Zf 

l/v2 

0 

0 

-1/  /2 

0 

i//2 

-1/  /2 

0 
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and  the  equation  (8)  can  be  reduced  to 


XI  =  (ci]^  +  d^i) 

2 

X4  =  Xi  (cii  -  dll) 
2 

yi  =  X;i  (C21  +  d2i) 
2 

74  =  '^1  {^21  -  <321  ) 

2 
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3.0 


FERRITE  MATERIAL 


Choice  of  a  ferrite  material  for  our  application  involves 
the  following  considerations. 


(a)  The  material  should  be  reproducible  on  a  piece  to 

piece  basis  -  particularly  free  from  density 
variations  and  4irMg  variations.  • 

(b)  The  material  should,  when  used  in  the  device, 
have  acceptably  low  rf  magnetic  loss  at  both  low 
signal  strength  and  high  signal  strength  (i.e., 
not  exhibit  nonlinear  loss  mechanisms). 

(c)  The  material  should  have  a  low  temperature 
coefficient  of  magnetization  -  this  normally 
implies  a  Curie  temperature  in  excess  of  300OC. 

(d)  The  ratio  Y4iTMs/t'^  should  be  between  0.4  and  0.8. 
Below  0.4  the  phase  shift  per  unit  length  becomes 
too  low,  above  0.8  the  ferrite  exhibits  low  field 
loss. 


Of  the  basic  families  of  ferrite  materials  that  are 
available  (hybrid  garnets,  magnesium  manganese,  spinels, 
lithium  ferrites  and  nickel  ferrites)  nickel  ferrites  are 
most  useful  at  the  frequency  range  and  power  level  in 
consideration.  They  have  high  Curie  temperature  (greater 
than  450°C),  suitable  magnetization  (1500  Gauss  to  3100 
Gauss)  and  acceptable  magnetic;  loss.  C-20  nickel  ferrite 
material  from  Countis  Industries  was  used  for  circulator 
design . 
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4.0  PHASE  SHIFTER 


4 . 1  Theoretical  Analysis 

Theoretical  analysis  was  performed  on  dielectric  and 
ferrite  loaded  ridge  waveguide.  The  geometry  under  study  is 
shown  in  Figure  3.  Dielectric  rib  is  fully  extending  across 
the  ridge  gap,  and  ferrite  slabs  are  partially  extending 
across  the  ridge  gap.  Exact  solution  to  the  given  problem 
is  very  difficult  to  obtain,  if  not  impossible.  To  approximate 
propagation  constant  of  the  dominant  mode,  transverse  resonace 
equation  was  derived  employing  ABCD  matrix  method. ^ 


In  ABCD  matrix  method,  the  electric  field  in  the  z 
direction  and  magnetic  field  in  the  y  direction  at  x  =  0  are 
related  to  the  electric  field  in  the  z  direction  and  the 
magnetic  field  in  the  y  direction  at  x  =  a  by  an  ABCD  matrix  as 


'Hy  (X  =  a)' 

CQ 

< 

1 _ 

Ez  ( X  =  a  ) 

Q 

U 

_ » 

Hy  (X  =  0) 
Ez  (X  =  0) 


(1) 


Because  the  electric  fields  at  the  waveguide  walls 
are  zero,  equation  (1)  can  be  written  as: 


Hy  (X  =  a) 

”  A 

B 

"Hy  (X  =  0) 

0 

_C 

D_ 

_  0 

l"TE-Mode  Solutions  for  Partially  Ferrite  Filled  Rectangular 
Waveguide  Using  ABCD  Matrices",  W.P.  Clark,  K.H.  Hering, 
D.A.  Charlton,  MTT  1966. 


12 


or, 


Hy  (x  =  a) 
0 


A  Hy  ( X  =  0 ) 
C  Hy  (X  =  0) 


(2a) 


For  nontrivial  solution,  the  matrix  term  C  must  be  zero. 
Therefore,  the  transverse  -^onance  equation  requires 
C  =  0 

By  applying  the  ABCD  method  to  each  region  in  Figure  3, 


matrix  quantity  C  can  be  calculated  as  follov/s: 

r  ir  TP  ir 

A  B  A^  Bj  A2  B2  Ar.  3n 

D  ^  C2  D2  Cr; 

ABCD  matrix  for  each  region  can  be  represented 

i)  Dielectric  Region  (including  air) 

* 

A  =  cos  Kj(  ^ 


B  =  -j 


C  =-j 


Sin  Kx  ^ 


Sin  Kx  ^ 


where 


D  =  COS  Kx 
Kx  =  e  -  6^ 

i  is  the  width  of  the  region 

Ferrite  region  (Ho  and  Mo  are  in  +z  direction) 
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Where  jB^  iz  capacitive  junction  susc:  tance.  The 
value  of  B  can  be  obtained  from  the  equation  LIow.^ 


r  /I  -  n  +  a 


(a  +  l/o.)/z 


1  -  a 


(A  +  A"  +  2C) 
AA'  -  c2 


jp  \  2  /I  -  o(\4rt  /5a  2  _  '  4  y2Q 


1  -  a2 


•1  -t-  1  +  /I  -  (b/  'g)2  1  + 

1  -  J  1  -  /l  -  (b/  \qy^  I  _ 


Where  A  = 


3  a2 


( iSa  ) 


1  +  a\2/a  1  +  ^l-(d/\q)  2 
1  -  a  j  1  -  /l-  (d/>.g)2 


3  +  .2 

1  -  -.2 


(15b) 


/4a  \2 

“  1  1  -  a2) 


(150 


(ISd) 


\g  = 


(15e) 


Yol  = 


(15f ) 


K0^£0  -  32 


(15g) 


Note;  The  value  of  Be  obtained  from  equation  (14)  assumes 

homogeneous  loading  of  dielectric  in  the  ridge  waveguide. 
In  case  of  inhomogeneous  loading  of  dielectric  in  the 
ridqe  wavea'.'ido,  the  value  of  Be  may  be  v'-ed  as  crude 
a  p  p  r  o  ;<  i  :r.  a  c  ion  . 


'"Wavenuide  H.indbook",  N .  Marcuvitz,  McGraw-Hill,  1951. 
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The  solution  for  the  propagation  constant  of  the 


geometry  shown  in  Figure  3  was  programmed  for  the  Hewlett 
Packard  9845  desk  top  computer.  A  copy  of  that  program 
is  included  as  Appendix  B  of  this  report. 

In  the  preceeding  analysis  it  was  assumed  that  the 
ferrite  region  was  fully  filled.  In  practice  it  will  only 
be  partially  filled  in  a  high  power  circulator.  It  was  found 
that  an  averaging  technique  could  be  used  to  adequately 
represent  the  partially  loaded  ferrite  structure.  The 
approximate  values  were  determined  by  a  set  of  iterations 
comparing  analysis  with  measured  results. 
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4 . 2  Phase  Tracking 

To  implement  the  design  of  a  four  port  differential 
phase  shift  circulator  using  two  magic  tees  and  two  non-identical 
nonreciprocal  phase  shifters  as  shown  in  Figure  1(d),  the 
following  conc.ltions  are  imposed  upon  the  phase  shifters: 

=  i^BF  ( 16  ) 

*^AR  "  ^BR  “  180°  (17) 

where:  is  insertion  phase  of  Phase  Shifter  A  in  forward 

direction. 

5z(ap  is  insertion  phase  of  Phase  Shifter  B  in  forward 
di rect ion. 

is  insertion  phase  of  Phase  Shifter  A  in  reverse 
direction. 

is  insertion  phase  of  Phase  Shifter  B  in  reverse 
di rect ion. 

Oue  to  the  frequency  dependence  of  phase  character¬ 
istics,  these  requirements  are  difficult  to  meet  over  wideband. 

There  are  many  parameters  which  can  effect  insertion 
phase,  such  as  width  of  the  dielectric  rib,  dielectric  constant, 
location  of  dielectric  rib,  dimension  and  location  of  ferrite 
slab,  and  biasing  magnetic  field  strength.  To  meet  the 
phase  tracking  requirement,  tiie  insertion  phase  differential 
of  the  two  phase  shifters  at  demagnetized  state  should  be 
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1 

I 

j 

reasonably  uniform  over  frequency  range,  preferably  about 
90°.  The  choice  of  the  dielectric  constant  and  the  width  of 
the  dielectric  rib  have  the  most  significant  effect  in 
setting  the  proper  amount  of  insertion  phase  differential. 

Typical  insertion  phase  characteristics  of  dielectric  and 
ferrite  loaded  phase  shifter  in  WRD750D24  double  ridge 
waveguide  are  shown  in  Figures  5  and  6.  Two  phase  states 
due  to  oppositely  directed  biasing  magnetic  field  are  shown 
along  v/ith  demagnetized  phase  state.  Here  we  define  longer 
phase  state  as  6  "  state  and  shorter  phase  state  as  3  state. 

The  insertion  phase  characteristics  of  an  empty  WRD750D24  double 
ridge  waveguide  is  also  shown  in  the  figures.  In  Figure  5, 
two  phase  shifters  having  the  same  dielectric  constant  are 
compared  with  respect  to  the  variations  in  the  width  of  the 
dielectric  rib.  In  Figure  6,  two  phase  shifters  having  the 
same  width  of  dielectric  rib  are  compared  with  respect  to 
the  variations  in  dielectric  constant. 

The  phase  tracking  requirement  given  in  equations 
(16)  and  (17)  can  be  satisfied  if  we  design  the  two  phase 
shifters  in  such  a  way  that: 

i)  The  insertion  phase  of  one  phase  shifter,  say 
phase  shifter  h,  in  B"*"  state  is  equal  to  the 
insertion  phase  of  another  phase  shifter, 
say  phase  shifter  B,  in3~  state  over  the 
operating  frequency  range. 

ii)  The  insertion  phase  of  phase  shifter  A  in 

•”  state  is  180°  degree  longer  than  that  of 
phase  shifter  B  in  .d"  state  over  the  operating 
frequency  range. 
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Dimensions : 

Wide  Rib :  = 

G3  = 

Narrow  Rib:Gi  = 
G3  = 


.014, 

Wd  = 

.035, 

G2  = 

.055, 
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Figure  5.  Phase  Characteristics  With  Variations  in  the  Width  of 
Dielectric  Rib  (Dielectric  Constant  =  13) 
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Phase  Characteristics  With  Variations  in  the  Dielectric 
Constant  (Width  of  Dielectric  Rib  Wn  =  .027") 


It  can  be  accomplished  by  employing  one  phase 


shifter  with  shorter  but  wide  dielectric  rib  (or  higher 
dielectric  constant)  and  another  phase  shifter  with  longer 
but  narrow  dielectric  rib  (or  lower  dielectric  constant). 

The  difference  in  length  between  these  two  sections  is 
made  up  by  length  difference  in  the  waveguide  employed. 

The  effect  due  to  the  insertion  of  matching  transformers 
can  be  minimized  by  designing  transformers  having 
approximately  the  same  insertion  phase  length  for  both 
phase  shifters.  Small  insertion  phase  difference  due  to 
the  matching  transformer  can  be  easily  corrected  by  adjusting 
other  parameter,  such  as  location  of  dielectric  rib  and  biasing 
magnetic  field  strength.  Neglecting  matcliing  transformer, 
actual  insertion  phase  differential  of  the  two  phase  shifters 
at  demagnetized  state  becomes: 


where 


is  the  phase  length  of 
at  demagnetized  state 

0g  is  the  phase  length  of 
at  demagnetized  state 


is 

the 

phase  length  of 

Lvj 

i  s 

the 

physical 

length 

hN 

is 

the 

physical 

length 

La 

is 

the 

physical 

length 

Ln 

= 

Lv;  + 

La 

wide  phase  shifter  per  inch 

narrow  phase  shifter  per  inch 

empty  normal  waveguide  per  inch 
of  wide  phase  shifter 
of  narrow  phase  shifter 
of  empty  normal  waveguide 


* 
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By  making  adjustment  of  length,  optimum  result  for  uniform 
insertion  phase  differential  can  be  obtained. 

Once  insertion  phase  differential  is  set  within  the 
acceptable  range,  fine  adjustment  can  be  made  by  adjusting 
location  of  dielectric  rib  and  ferrite  slab,  and  biasing 
field  strength.  Typical  differential  phase  characteristics 
between  “  state  and  ?  +  state  with  respect  to  the  location 
of  the  dielectric  rib  is  shown  in  Figure  7(a)  and  (b). 

When  the  dielectric  rib  is  placed  near  the  center  of  the 
ridge,  it  gives  more  differential  phase  at  higher  frequency. 

When  the  dielectric  rib  is  moved  toward  the  side,  it  gives 
more  differential  phase  at  lower  frequency. 

In  practice,  breadboard  phase  shifters  were 
built  based  on  computation,  and  absolute  phase  measurement 
was  performed  with  variations  in  the  width  of  dielectric 
rib,  dielectric  constant,  location  of  dielectric  and  ferrite, 
and  Hq(3  field  strength.  A  simple  computer  program  was 
written,  and  measured  phase  data  were  fed  into  the  computer 
program  to  find  optimum  geometry.  Introduction  of  matching 
sections  introduced  small  deviation  from  optimum  result,  but 
compensation  was  done  by  adjusting  the  location  of  the  dielectric 
rib . 
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Figure  7 (a) . 


Differential  Phase;  Dielectric  Rib  is  Loaded  at 
the  Center  of  the  Ridge  (WRD750D24) 


Gl 

G3 


.069,  Wp 
.019,  Tf 


.035,  G2  =  0,  Wf  =  .050 
.030,  =  13 


J - 1 _ I _ I _ L 


Figure  7(b).  Differential  Phase;  Dielectric  Rib  is  Loaded  Near 
the  Edge  of  the  Ridge  (WRD750D24) 

=  .014,  Wp  =  .035,  G2  =  .055,  Wp  =  .050 

Gj  =  .019,  Tf  =  .0-30,  = 


GH 


13 


Figure  3  shov/s  the  actual  dimension  for  phase  shifters 


and  experimental  results  on  phase  are  shown  in  Figure  9. 

(Note:  Sign  is  inverted  on  phase  scale,  that  is,  bottom 

side  has  longer  phase). 

4 . 3  Matching 

Chebyshev  type  3  section  dielectric  loaded  quarter 
wave  transformers  are  utilized  for  both  phase  shifters.  As 
mentioned  above,  the  insertion  of  the  matching  transformers 
introduced  small  deviation  from  optimum  result.  By 
adjusting  the  location  of  the  dielectric  rib  slightly, 
optimum  result  was  obtained.  Results  on  VSV7R  for  both  phase 
shifters  are  shown  in  Figures  10  and  11. 

4 . 4  Insertion  Loss 

Figures  12  and  13  show  total  insertion  loss  for  both 
phase  shifters.  These  measurements  include  the  line  lengths  of 
waveguide  employed.  Insertion  loss  can  be  reduced  to  some 
degree  by  improving  VSWR  and  using  low  loss  bonding  material. 


PHASE  SHIFTER  A 


PHASE  SHIFTER  B 


.010" 

Wq  (Width  of  Dielectric  Rib)  .035" 
G2  .058" 
Wp  (Width  of  Ferrite  Slab)  .050" 
G3  .020" 
Tp  (Thickness  of  Ferrite  Slab)  .030" 
Length  3.890" 
Dielectric  Constant  13 


.ij25" 

.0225" 

.0455" 

.050" 


.030" 


.030" 


5.090 

16 


Figure  8.  Dimensions  for  Phase  Shifters 
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5.0 


NIAGIC  TEE 


5 • 1  General  Considerations 

The  major  effort  associated  with  the  design  of  the 
hybrid  magic  tee  is  the  impedance  matching  at  the  junction 
over  the  intended  bandwidth.  Methods  of  matching  the 
junction  include  the  use  of  reactive  elements  at  the  junction 
cavity,  and  transformers  between  the  flanges  of  each  arm  and 
the  actual  junction  at  the  cavity  wall  for  each  arm.  In 
addition  to  matching  problem,  there  are  other  considerations 
such  as  dissipated  losses,  phase  balance,  and  pov;er  handling 
capacity.  To  minimize  dissipated  losses,  it  is  desirable 
to  avoid  the  use  of  dielectric  materials,  bonding  agents 
or  screws  inside  the  waveguide  cavity.  Phase  balance  can 
be  obtained  as  long  as  physical  symmetry  is  maintained.  In 
terms  of  pov/er  handling  capacity,  it  is  desirable  to  avoid 
sharp  capacitive  obstacles  or  low  height  construction  in 
the  waveguide. 

5 . 2  Computer  Analysis 

In  an  attempt  to  study  and  design  magic 
tee  using  double  ridge  waveguide,  a  computer  program  was 
written  to  find  cutoff  frequencies  and  characteristic 
impedance  of  double  ridge  waveguide.  Cross-sectional  shape 
of  doufc  le  ridge  waveguide  and  its  equivalent  circuit 
representation  is  shown  in  Figure  14  and  equations  used 
for  computation  are  summarized  below  without  derivation. 
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a 


4.  Do'..jble  RiJae  Wavequide  Cross-Section  and  Its 
iiiuivalont  Circuit  Representation 
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(Original  derivation  of  these  equations  can  be  found  in 
Reference  3 ,  4 ,  5 . ) 

The  equations  which  govern  the  cutoff  conditions 
of  the  TE^q  type  of  nodes  are  given  by: 
b 

cotk^l  -  -  tan  s/2  -  B/Yq]^  =0  (1) 

d 

b 

cotk^l  +  -  cotk^  s/2  -  B/Yq]^  =0  (2) 

Equation  (1)  applies  to  the  odd  TE^q  nodes  and  (2)  applies 
to  the  Oven  TCp,.^  nodes.  is  the  propagation  constant 

in  the  x  direction  at  the  cutoff  and  is  liven  by  k^  =  I'/'-q. 
The  characteristic  admittances  in  the  transverse  direction 
( x-di rect ion )  Y^|  and  Yq2  are  defined  as: 


Yol 


1 

b 


1 

d 


(3)  S.B.  Cohn/  "Properties  of  Ridge  Waveguide",  Proc.  IRE, 

Vol,  35,  pp.  783-788;  August,  1947. 

(4)  T.C.  Mihran,  "Closed  and  Open-Ridge  Waveguide",  Proc.  IRE, 
Vol.  37,  pp.  640-644;  June,  1949. 

(5)  S.  Hopfer,  "The  Design  of  Ridge  Waveguides",  IRE.  Trans., 
Vol.  MTT-3,  pp.  20-29;  October,  1955 


The  value  of  the  normalized  susceptance  term  B/YqI,  which 
represents  the  effect  of  the  step  discontinuity,  can  be 
obtained  from  IVaveguide  Handbook,^. 

Employing  power-voltage  impedance  definition, 
the  characteristic  impedance  at  infinite  frequency  for 
the  TEj^g  mode  is  given  by: 


2  *  120:t2 


^pv 


T  “  2 

2"  —  cos  +  ' 

e  2  c 

o 


2  sin2  2 


2cos 


2  . 


2d 


b*  s in' 


sin2 

4 


where 


'q  is  cutoff  wavelengtii  of  ridge  waveguide. 
1  =  ■■  (a-s) 


-  o 

—  —  —  '^'Oc;h 

1+ ( d/b ) ^ 

4(d/b) 

V..  V../0  1 1 

(-1 

_l-(d/b)2_ 

———————— 

l-(d/b)2 

'b' 

_ 
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The  characteristic  impedance  at  any  frequency  f  is  obtained 


5 • 3  Magic  Tee  Design 

Initially,  magic  tee  was  designed  based  on 
intuitive  reasoning.  At  the  cavity  junction  impedance  ratio 
of  H-arm,  colinear  arm,  and  E-arm  is  assumed  to  be  1:2:4.  Then 
Chebyshev  type  multisection  quarterwave  transformer  between 
flanges  of  each  arm  (V?RD750D24)  and  cavity  junction  of  each 
arm  were  computed  using  equation  (4).  Cutoff  frequencies  of 
TE]_g  and  TE20  node  at  the  junction  were  also  computed.  Based 
on  the  computed  results,  two  sets  of  brassboard  magic  tees 
were  made.  Brassboard  was  chosen  so  that  additional  metal 
could  be  soldered  in  place  v/henever  design  modification 
is  necessary.  Figure  15  shows  simplified  sketch  of  the  magic 
tee.  For  the  purpose  of  evaluation  of  characteristic  impedance 
of  double  ridge  waveguide  and  verification  of  the  theoretical 
result,  each  magic  tee  was  made  in  three  separate  parts. 

Part  1:  H-Arm 

Designed  with  I'^ur  section  step-down  quarterwave 
transformer  beuv'een  WRD750D24  and  H-arm  junction. 
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Part  2: 


E-Arm 


Designed  with  3  section  step-up  quarterwave  transformer 
between  WRD750D24  and  E-arm  junction. 

Part  3:  Colinear  Arm  Including  Junction  Cavity 

Designed  with  4  section  step-down  quarterwave  transformer 
between  WRD750D24  and  colinear  arm  junction. 

Since  Part  1  and  Part  2  of  the  magic  tee  do  not  include  junction 
cavity,  these  parts  were  used  to  verify  computed  characteristic 
impedance  of  the  double  ridge  waveguide.  VSWR  of  H-arm  (two 
H-arm  in  back  to  back  configuration)  was  measured  and  Smith  Plot 
of  the  measured  VSWR  is  shown  in  Figure  16.  Smith  Plot  of  VSWR 
for  the  same  configuration  was  computed  using  transmission  line 
equation  employing  characteristic  impedance  given  by  Equation  (4). 
Computed  results,  shown  in  Figure  17,  were  in  good  agreement 
with  experimental  results. 

Due  to  the  difficulties  of  transmission  line 
model  including  junction  cavity,  no  attempt  has  been  made  to 
analyze  magic  tee  as  a  whole  theoretically.  Instead,  it 
was  intended  to  find  optimum  matching  structure  experimentally 
by  using  reactive  element  at  the  junction  cavity.  The 
effect  of  reactive  element  can  be  measured  at  the  flanges 
of  each  arm,  and  then  these  measured  impedance  can  be 
transformed  to  cavity  junction  wall  using  transmission 
line  equation.  Based  on  these  transformed  impedance,  necessary 
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modification  can  be  made.  This  method  requires  iterative 
process  of  microwave  measurement,  engineering  analysis,  and 
modification.  However,  this  method  has  not  been  pursued 
due  to  the  time  constraint.  Therefore,  magic  tees  designed 
by  Raytheon  for  previous  projects  were  used  for  final 
assembly  of  circulator.  These  magic  tees  were  designed 
for  frequency  range  between  8  to  18  GHz.  Performance 
data  of  magic  tees  are  shown  in  Figures  18  through  25. 
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6.0 


HIGH  POWER 


High  power  measurement  has  not  been  performed  for  the 
circulator  assembled.  According  to  estimation  based  upon 
available  informations,  circulator  v^ould  require  pressurization 
and  efficient  cooling  for  safe  operation  at  the  specified 
high  pov/er  level. 

6 . 1  Peak  Power 

6.1.1  Magic  Tee 

Magic  tee  used  for  circulator  assembly 
include  capacitive  obstacles  for  optimum  matching.  In  order 
to  guarantee  safe  operation  at  the  specified  power  level, 
pressurization  of  waveguide  would  be  necessary.  Pressurization 
windows  for  the  specified  power  level  do  not  exist  at  this  time. 
Microwave  Research  Corporation  Model  270  is  rated  at  10  KW 
peak/600  watts  average  for  a  VSWR  of  <1.1.1  and  loss  of  <0.1  dB. 
Microwave  Development  Labs  have  a  unit  under  development.  Further 
development  will  be  necessary  if  the  desired  power  levels  of  20  KVv’ 
peak/1  KW  average  are  to  be  met. 

6.1.2  Phase  Shifter 

There  are  two  effects  of  concern  in  the 
phase  shifter  section,  breakdown  and  nonlinearity.  Nonlinearity 
is  the  increase  in  insertion  loss  above  a  certain  peak  power 
threshold.  Ferrite  material  used  for  phase  shifters  design 
has  been  utilized  previously  for  the  similar  frequency  range 
with  comparable  geometry.  To  prevent  peak  power  breakdown, 
air  gap  should  be  removed  in  bonding  of  dielectric  rib. 
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6 . 2  Average  Power 


6.2.1  Magic  Tee 

The  average  power  capability  of  the 
magic  tees  would  be  adequate  if  conduction  cooling  to  liquid 
cooled  phase  shifters  is  used. 

6.2.2  Phase  Shifter 

Since  the  phase  tracking  is  essential 
for  the  operation  of  the  circulator,  efficient  cooling  of 
phase  shifter  is  required  to  prevent  phase  change  due  to 
temperature  variations.  Liquid  cooling  with  pressure  drop 
of  10  psi  would  be  adequate. 


7 . 0  CIRCULATOR  ASSEMBLY  AND  TEST  DATA 

Cross-sectional  view  of  phase  shifters  for  the  circulator 
assembly  is  shown  in  Figure  26.  Final  data  on  circulator 
assembly  are  shown  in  Appendix  A.  WRD750D24  calibration 
standards  were  used  in  the  calibration.  The  designation 
Port  X-Y  indicjates  that  Port  X  was  the  input  port  while 
output  signals  were  measured  at  Port  Y. 
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Figure  26.  Cross-Sectional  View  of  Phase  Shifters  for  ('irculator 
As.semljly  (Port  1  is  Toward  Viewer) 


8 . 0  CONCLUDING  REMARKS 

Four  port  cli f ferential  phase  shift  circulator  under 
Contract  No.  N00173-80-C-0434  for  Naval  Research  Laboratory 
has.  been  built.  For  frequency  range  between  8  to  16  GHz, 

11  dB  minimum  isolation,  1.8  dB  maximum  insertion  loss,  and 
maximum  VSWR  of  2:1  were  achieved.  Even  though  the  performance 
of  the  circulator  designed  did  not  meet  the  originally 
intended  specification,  it  has  demonstrated,  at  least,  the 
feasibility  of  building  the  broadband  four  port  differential 
phase  shift  circulator. 
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APPENDIX  A 
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appendix  b 

COMPQTER  PROGRAM 


B1 


▼ 


1 0 

PEM  Pro'gram 

"PIEGE" 

'  ; 

20 

OPTION  ERSE 

1 

:'U 

EIH  Nbt'  20:' 

,  Pb»  '  20 

' ,  Obt  ■;  20  > ,  Pbo 20  > 

4  0 

CON  R'7',E': 

7  ■  ,  F  ,  K  ,  1 

k  ,  Gama,  Sat  ,  Rem,  E,  E 

50 

C  r  1  t  = . 000 1 

60 

Ep£= . 000 1 

7  0 

Ea=l 

'z‘  0 

M  =  0 

90 

N  =  0 

1  0  0 

Rl=. 253 

1  10 

R2= . 0 1 4 

II'M  m.::=.£i35 
130  R5=.050 

140  He:=.0;'0 

150  R4=  .  173-'  R3  +  R3  +  R5  +  R(;  ■ 

li;0  R7=.25'? 

170  E=.321 

130  D=.l375 

130  El =.030 
200  E2=D-2*E1 

210  Ef=11.4 

220  El=l 
230  E2=l 
240  E3=13*.7S 

250  E4=l 

2';0  E5='2^E1+E2'  '.2fEl  Ef+E2  Ea' 

27  0  E  €■  =  1 

2S0  E7=l 
230  iJ(:'n'i  =  0 
300  Gam a =3. 5 

310  Sat =1320 
320  P€mf;'=l320 

3  3  0  R 1=  m  r'l  =  -  P  •;  rn  p 

340  Fi=3 
350  Ft =16 
360  Ds1f=l 

3  70  R  1  '  =  R  1 

330  R':2:'=R2 

330  R'3'=R3 

400  R  4  '  =R4 

4  10  R':5'=R5 

420  R'-6:'=R6 

4  3  0  R  7  .'  =  R  7 
440  E'  1  '=E1 
450  E':2.'=E2 

460  E '  3  '  =  E  3 
470  E'.4;'=E4 

430  E'5:>=E5 

4  3  y  E  '  6  '  =  E  6 

500  E'::7:'=E7 

510  I  URGE  4X,  "R1  "  ,  7X,  "R2"  ,  7X,  "R3"  ,  7X,  "R4"  ,  7X.  "R5''  ,  7>:,  "R6"  .  7N:.  ''R7" 

520  PRINT  USING  510 

5  30  IMRGE  2:',  7'- DD.  DEED.  2X  ' 

540  PRINT  USING  5 30 ; R 1 , R2 , R 3 . R4 . R5 , R6 , R7 

550  IMRGE  4>:,  "El"  .  7X.  "E2"  .  77,  "E3"  ,  77,  "E4"  .  7;;.  "E5"  .  77.  "E6"  ,  7:;,  “E7" 

560  PRINT  USING  550 

570  PRINT  USING  5 30 ; E 1 , E2 , E 3 , E4 . E5 , E6 , E7 

5  3  0  I  M  R  ij  E  4 . . ,  "  E  "  ,  cl  ,  "  El "  ,  3 .  ^  "  Ei  1  "  ,  r  ,  "  E  t  "  ,  6  -s ,  "  u  am  a ",  6  .s ,  "  3*  at  ",  6  3! ,  "  P  €  m  +  " 

530  PRINT  USING  530 

600  IMRGE  27, 5'  EE.  EEEE,  2::  ' ,  7, 2' EEEE.  E.  37  > , 

6  1  0  PR  I  NT  US  I  NG  600  ;  E  ,  E  ,  E  1  ,  Et  ,  Gama,  S  at. ,  Pemp 

m 


F  =  Fi 

b  J  0  R  S'  n  I  =  0 

640  L  b 1 :  R  b  t  =  6 . S 

6  50  I b  t  =  0 

660  IF  Sat =0  THEN  GOTO  Lb2 

670  IMfiGE  ,  "RerM="  ,  MDDDD,  5X,  ■•**«*#*#♦**■•,  , 

6S0  PRINT  USING  670; Rem 

6  9  0  R’  e  m  =  R  e  m  ,  5  *  t  1  +  E  2  E  2  +  2  E  1  > 

7  £t  £'  L  b  2  ;  K  =  2  ■»  P  I  *  F  ■  1  1  .  S 
710  Kk=K*K 

7  2  0  CALL  N  e  w  X  o  n  R  b  t  ,  I  b  t.  ,  C  r  i  t  ,  E  fj  s  > 

730  M=N+1 

740  IF  Rem=0  THEN  Lb3 

750  GOTO  Lb4 

760  Lb3:  RbO'::M;'=Rbt 
770  Bet a=Rbt 
7S0  GOTO  Lb5 

7  9 1'  L  b  4  :  D  b  t.  =  R  b  t  -  R  b  o  ■  M  >  :>  -fr  2  *  B  1  D 

8  y  0  E  €  t  .=1= R!  b  o  *■.  t'1  >  +  D  b  t 

310  Lb5:  Phase=Beta*lS0  PI 

820  IMAGE  DD.DD,''  GHz  "  ,  5;': ,  K  Y=  "  ,  D .  DDDDE  ,  "  I  NCH  "  ,  5X  .  "  PHASE=  ".D.DDDDE."  DEG  I 
NCH" 

8;  3  0  PRINT  U  S  I  N  G  3  2  O  ;  F  .  E  e  t  a ,  P  h  a  =  e 

340  IF  Rern  =  0  THEN  Lb7 

3  5  U  IF  R  e  rn  ■  0  THE  N  L  b  6 

860  Pbr  (ri  ::'=Phase 

870  GOTO  Lb 8 

830  Lb6:  Nbt  ■..M>=  Phase 

390  GOTO  LbS 

900  Lb7:  Obt  ■:N  ,)=Phai  e 

910  Lb8:  IF  F=Ff  THEN  Lb9 

920  F=F+Delf 

9  3  0  G  0  T  0  L  b  2 

940  Lb9:  IF  Rem  0  THEN  Lb  10 

9  5  0  L  g  =  7  2  0  ‘  0  b  t  1  j  +  0  b  t  ■.  M  ■  ■  4 

9 6  0  I N  AGE  " 0  u  ar  t  e  r  W  a" e  L  e  n g  t  h  = " , D . D  D  DDE."  I  N  L  H  " 

9  7  0  F'  R'  I  N  T  U  S’  I  N  G  9  6  O  ;  L  g 

930  Lb  10:  IF  Sat=0  THEN  Lb99 

990  N=N+1 

1000  ON  N  GOTO  Lbl 1 , Lbl2, Lbl4 

10  10  L  b  1  1  :  R'  e  m  =  R  e  m  fi 
1020  GOTO  Lb  13 

10  30  L  b  1  2  :  R'  e  m  =  R' s'  rn  r 
1  0  4  Cl  L  b  1  3  :  F  =  F  1 
1  0  5  0  M  =  0 
1060  GOTO  Lbl 
1070  Lb  14:  M=1 
1030  F=Fi 

1090  PRINT  LIN'  1  "DIFFERENTIAL  PHASE'- DEG  INCH'" 

1100  PRINT  "FREQ'  GHz:'  E-  -  Bo  E+  -  Bo  E-  -  E+" 

1110  L  b  1  5 :  D  n  R  =  N  b  r  ■ :  M  '  -  P  b  t  ■:  M  ,■ 

1  1  20  Drio  =  Nbt  '  M  '  -Obt  '  M  ' 

1130  Dp o  =  P b  t  ':  H  - 0 b t  M  ) 

1  1  40  IMAGE  K,  DD.DD,  7:";,  MDDD.  DD,  5::,MDDD.  DD  .  57; ,  MDDD  .  DD 

1  1  50  PR  I  NT  US  I  NG  1  1  40  ;  F  .  Dno ,  Dpo  ,  Dr.p 

1160  M=M+1 

1170  F  =  F  +  D  e  1  t 

1130  IF  F::=Ff  THEN  Lb  15 

1190  Lb99:  STOP 

1200  END 

1210  S  LI  E  .T  u  n  i;  t  '  E .  k  ■:  1  ,  E  ,  D  ■ 

1220  Ar =D  E 
12  30  L':=2*PI  K:1 

1240  Cb=  '  4*Ar  ■ 1 -Ar  2  '  '  2 
1250  Apb='.  1-'  D  L':  '  2  '  .  5 
1260  Apb='l+Apb'  'l-Apb' 

RT 


12  70 

Apb 

-Apb*'  1  +  A  r  ' 

1  - 

Ar  ;•  ..' 

■  2  Ar  ' 

1  230 

Apb 

=  Apb+  '.  3  +  Ar  ■  2 

n  U 

-  Ar  ■ 

•  ^  .1 

1290 

Ab  = 

';  1-';  B  L'Z  '  -  2  ' 

Cj 

1  300 

Ab  = 

1  +Ab  '  1  -Ab 

1  3  1  0 

Ab  = 

Ab*  ':  ',  1  +Ar  ,'  ■  < 

1-Ar 

>  .  -v 

2*Ar  .' 

1  320 

Ab  = 

Ab-  1  +3*Ar  -2 

>  ^  1 

-Ar 

•  2^  j 

1  3  30 

Bl  = 

1  -A(  -2  <  4*  Ar  > 

1  340 

Bl  = 

B  1  ♦  V  1  +Ar  ) '; 

1  -Ar 

}  '•  ■ 

'.  Ar+l  Ar  .'  2> 

1  550 

Bl  = 

LOG'.  El  :> 

1  360 

B2  = 

2+  Ab  +  Apb  +  2* 

Cb  3 

Hb* 

Apb-Cb 

2  -> 

13  70 

B3  = 

( ’  ■■  'I  4*Lc  >  '  '  2 

1  330 

B3  = 

F  .*'..(  l-Ar).- '. 

1  +Ar 

)  .1  •  ••  •' 

4  ♦  A  r  ' 

1  390 

B3  = 

3  *  ':  (  5  *  A  r  ■  ■  2  - 

1  ) .  '■ 

1  -Ar 

■■■■2  +  4*C 

b*Ar  2 

1400 

B.:  = 

*E..  L'l*'  El+E 

2  +  B3 

> 

14  10 

SUI 

.HD 

1420 

3  U 

H  s  1...I  r.  fi ':  R  b  t  , 

Ibf  , 

C  r  1  f 

,Eps:' 

14  30 

1  =  1 

1440 

Id  = 

1 

1450 

Lnl  : 

R3-  =  Rbt 

1  460 

I  3  = 

Ibr 

1  4  70 

IF 

It't.  =0  THEN 

Ln2 

1  4  3  0 

I  3  = 

-Ibr 

1  4  90 

Ln2 : 

R.:'  =  R3 

1  5  0  0 

Io  = 

I  3 

1510 

Ln3: 

ON  Id  GOTO 

Ln4 

.  Ln5 

I  520 

L  r"'  4  ; 

Pds 1 =Ep3 

15  30 

Ids 

1  =0 

1540 

GOTO  Ln6 

1550 

Ln5: 

Pds I =0 

1560 

Ids 

1  =  Ep3 

15  70 

Id  = 

1 

1530 

Ln6: 

Rr"  =  R':'-Pds  1 

1590 

Ir"  = 

I'O-Ids  1 

1  6  0  0 

Pp  = 

P'j  +  Rds  1 

1 6  1  0 

Ip  = 

I'j+Ids  \ 

1620 

Pbt 

=  P.:. 

1  6  30 

I  bt 

=  Io 

1 1:'  4  0  CALL  F  ij  r'l  d  '  P  b »  .  I  b  f  ,  F'  r  ,  I  r  ■ 

1650  C  M  L  L  L  ab  i  '  P  r  ,  Ir  ,  C,  r  • 

1660  IMPGE  MD.  DDDriE,  4,':.  MD,  DIiDDE.  6::.  D.  DIiriliE 

16Si0  IF  Cr‘  Cf  1  T'  THEN  R e- 1  _fti a i  fi 

1  6 '?  0  P  t  o  =  F'  f 

1700  Ifo=Ir 

17  10  R  b  t  =  R  m 

1720  I  b  t  =  I  rn 

1  730  CALL  Fund'.  Pbt  ,  Ib*  ,  Pn,  I  r  :■ 

1740  R  f  re  =  P  r 
1750  If  rn  =  I  n 
1760  R  b  t  =  R  p 
1770  I b  T  = I p 

17Sti  CALL  Fund'  Rbt.  ,  Ibf  ,  Pr,  Ir 

1  7 '?  0  P  f  p  =  R  r 

1300  Ifp=Ir 

1  '3  1  01  R  f 'd  =  R  f  o  “  R  f  rit 

1320  I  f  d=  1  f 'i- I  f  m 

1  3  30  CALL  Cd  i  "  ■:  Pf  d  ,  I  f  d  ,  Rde  1  .  I  de  1  ,  Rf  d  ,  I  f  d  ) 

1  3  4  £t  R  f  =  =  R  f  m  +  P  f  fj  -  2-^  R  f  i:i 

1  3  50  I  f  a  =  I  f  f"  +  I  f  p  -  2  ■*  I  f  III 

1  360  CALL  Crnu  1  t  ■  Pde  1  ,  I  de  1  ,  Pde  1  ,  I  d6  1  ,  R'd€-  Is,  Id 

1  370  CALL  C  d  i  ■' Pf  3  ,  I  f  s  ,  Rde  1  s  ,  I  ds  1  s  ,  Rf  s  ,  I  f  s  ■ 

1  3 3 0  C  ALL  C  abs  ■:  R f  s  ,  I  f  s  ,  C  f  s  :• 

1390  IF  Cf 3  .0000001  THEN  Ln9 
1900  OH  Id  GOTO  Ln7,Ln3 
1910  Ln7:  Id=2 

1920  IF  r  10  THEN  Lnll 
1950  1=1+1 


^Ab 


B4 


•.n  1.0 


1940 

ij  U  T  LI  L  fi  o 

1  950 

LfiS:  GOTO  Lnl0 

1  9  E  01 

Ln9:  CRLL  Crnult 

■:Rf.;., 

1970 

C  R  L  L  C  rti  i.j  1  t  •-  R  f  d  , 

Ifd,R 

1  930 

CRLL  Cdv.,":  Ros,  I 

,  Rd 

1  9  ■:*  0 

R  0  i  d  =  1  -  2  ♦  R  o  s  d 

2  0  0  0 

I.:.i-d  =  -2*Iosd 

2  0  1 13 

CALL  C  i  q  r  t  R  o  =  d 

,  Io£d 

2020 

R  •;  Cl  =  d  =  1  ~  F!  c  0  £  d 

2030 

I.:osd=-I.:c.id 

2  0  4  0 

CRLL  C  rn  u  1  r  '■  P  f  d  , 

I  f  d  .  R 

2  01 5  01 

CRLL  C  d  1  ' <  R  d  z  ,  I  d  z  ,  R » ' 

2  0  b  0 

Rdz=-Rdzz 

2070 

Idz=-Idzz 

2  0  S:  0 

CALL  Cabs< Rdz, Idz,Cd 

2  0  9  0 

IF  Cdz<. 0000001 

THEN 

2  1  00 

P  b  t  =  R  B  +  R  d  z 

2110 

Ibt.  =  I  =.+  Idz 

2120 

IF  I: 10  THEN  Ln 

1  1 

2  1  30 

1  =  1  +  1 

2  1  40 

GOTO  Lnl 

2150 

LnlO:  PRINT  "ZERO  HER 

2  1  30 

G  0  T  ij  R  B  t  rn  a  \  fi 

2  1  T'O 

Lnll:  IMAGE  "DI 

v'ERGE 

2130 

PRINT  U3ING  217 

0 ;  C  r 

2  1  9  0 

F'  €  t  fii  a  1  r'l  •  b  U  E  E  X  I  T 

2  2  0  0 

3 U BEND 

22  10 

iUB  C  F'lu  1  T  '  : :  1  ,  r'  1 

1  •  '<  4  1 

222  0 

::3  =  ::i+::2-ii*V2 

2  2  3  0 

'i'  3  =  i  1  ‘i‘  2  'i'  1 

2240 

SUBEND 

2  25  0 

SUB  Cdi'"::i,Yl, 

2  «  i’  2 

2230 

D  €'  r'l  0  fii  = ! !  2  2  +  i’  2 

r'  2 

2  2  7  0 

773=  1  7!  1  +  X2  +  V  1  ■»  Y2 

D€r. 

2230 

V  3  =  *  'r'  1  * !  !2  ~  .  .1  V2 

D«r, 

2290 

SUBEND 

2  3  0  0 

SUB'  C  ab i  '  X  .  V  ,  R m 

2  310 

F' fi'i  =  '  . ■*■,:.+  r  *  i  '  ,  5 

2  3  20 

SUBEND 

2  3  3  0 

SUB  C  =  1  n  '  1  , 'i' 1  . 

-'•i  2  <  'r  2 

2  340 

7:2  =  SIN':  :-:i  ■  +  ■  e:-:p 

Y  1  '  + 

2  350 

’r  2  ==  C  0  S  ]  Y 1  '  ^  ’  E  P 

■:  1  ■  - 

2  3  6  0 

SUBEND 

2  3  7  0 

'3  IJ  E  C-  c  C'  3  !  V 1  ,  '  i'  1  , 

-  Y  2  «  ’i*  2 

2  3  3  0 

:<2  =  C0S'  Xl  EXP 

Y  1  ■  + 

. D. DDDDE 


'i';  =  -SlN'  :;i  e::pvvi  :'-e:;p^  -'.'i  > 
SUEEHD 

S  U  B  C  s  q  r- 1  i  :  1  ,  Y  1  ,  :<  2.'t'2  :• 

H  ag  =  ■:  1  > : ;  1  +  Y  1  *  Y'  1  ,i  .25 

IF  THEN  Lbal 

IF  THEN  Lba2 

Plng  =  RTN  ■  'i'  1  1  2 

G  0  T  0  L  c  B  q 

Lb  al  ;  Rng=  (.  P  I  +hTN  ■  Y  1  : :  1  ■  ■  2 
GOTO  Lciq 

Lba2:  IF  ■i  l  0  THEN  Lba3 

Rrig  =  PI  4 
GOTO  Lcjq 
Lba3:  Hnq=-PI  4 


5  30 

Lcsq 

:  ::2  =  Mag  +  C0S' 

R  n  g 

540 

V'  2  — 

M  ag  +  S  I  N '  Rfig  j 

550 

SUBEND 

5  b  0 

SUB 

F  u  r'l  d  '■  R  b  1-  ,  I  b  t 

.  Rr- 

570 

DIM 

Rk  1  •  7  1 , Ik  :  ■  7 

■  .  R 

530 

D  I  M 

R'  fn  a  '■  11  '  ,  I  rn  a ' 

1  1  ' 

590 

COM 

A  '  +  '  ,  E  '  +  '  ,  F  , 

K  .  K 

1 


2800 

P0  =  0 

28  10 

11  =  1 

2  8  2  0 

CRLL  C  niu  )  f  '  Rbt  ,  I  bt  ,  Rbt 

,  I  b  t  ,  R  b  t 

3  ,  I  b  t  i 

2  8  5  0 

FOR  1=1  TO  7 

284  0 

Rk  :=kk»E' I'-Rbt 3 

2  8  5  0 

Ik,:=-Ibr3 

2880 

CRLL  C  3  q  r  t  <  P  k  x ,  I k  . . , R  k  x 

•  I  >  ,  I  k  X 

I  >  > 

28  7  0 

IF  lk:.x  l:':0  THEN  CRLL 

Neq^  Rkx( 

I  >  ,  I  k  X  I 

2  8  8  0 

NEXT  I 

289  0 

U=1 

2700 

Rxxv=0 

2710 

I xxy=0 

2720 

IF  S.at=0  THEN  Lbbl 

2730 

1.1=  1  o  +  1  “  '■  b  =tni it-  ■  V 

F  *  1  000  :■  :■ 

••  ..  Cj 

2740 

U  =  If  +  ' .  1  ~  LI  R  E  S  V  f'  1='  rn  ’  S  -it  > 

1 . 5 

2750 

IJ=1 

2780 

I  ;<  =  Ij  ifi  l  i'*  R'  €  Rl  1  F  *  1  U  U  y 

2770  Lbbi:  Rho  =  U  ( I  xxM  -  2  ' 

27SO  I  r  «  r  =1=  I  xxy  .  U 

2  7  9  £'  R  k  k;  =  K  k  *  E  (  5  >  R  ki  o  -  R  b  t  3- 

2  S  0  0  C  R  L  L  C  3  q  r  t  <.  R  k  x  ,  I  k  x  ,  R  k  x  ■;  5  ■  ,  I  k  x. 5  '  > 

28  10  IF  Ik.:(5>  0  THEN  CRLL  Neq  ■  Rk  x  5  '  ,  I  k  x  ■  75  ,  Pk  x  ■  ,  I  k  x  ■.  S  ..  . 

2820  Kxl=Rk::<i::' 

2  8  8  0  E  c  =  0 

2840  IF  Kxl<=0  THEN  Lbb2 
2850  IF  D=E  THEN  Lbb2 
2  8  8  0  CRLL  J  i.j  n  c  t  <  E  ,  k  x  1  ,  E  ,  D  > 

28  70  E  c  =  “ E  c  *  K  X  1 .  K 
2880  Lbb2;  FOR  1=1  TO  7 
2  8  9  0  R  k  X  i=  R  k  x  x  I  >  •«■  R  <  I  > 

2900  Ik  xa=Ik  ■<(  I  ::"»R':  I  ■ 


2  9  1  0 

CRLL 

c 

5  1  r*i  x  R  k  X  a ,  Ik  ■■■  a  *  P!  i  n ,  I  i  n  ■ 

2920 

R  i  n  X 

I  ‘ 

=  P:sfi 

29  30 

I  i  n  x 

I  ;■ 

=  Un 

2  9  4  0 

CRLL 

c 

'I  0  i  '■  P'  k  .  -  a ,  I  k  y:  a ,  R  *1  r'l ,  I  c  r*i  ■ 

2950 

Pcnx 

I 

=  R.;r-i 

2  9  8  0 

I n 

I  ■ 

=  I n 

2  9  70 

Pirn  a '■ 

I  • 

=  Px:r-, 

2  9  8  0 

I  Ti  l  a ' 

I  ' 

=  I  c  n 

2  9  9  0 

CRLL 

X 

rn  i.j  1  t  1  P  k  1  I  '  ,  I  k  '  I  '  ,  P  i  n ,  I  i  rc ,  P  k  i  n  ,  n  .  = 

3  0  0  0 

CRLL 

c 

rn  i.j  1  r  '  P  0  ,  I  1  ,  P:  k  y  =  ki ,  I  k  -  -  fi ,  P  III  b  ,  I  rn  b  ' 

3010 

P!  ni  b  '■ 

I 

=  -Rrnb  k 

30  30 

I  mb ' 

I  ' 

=  -  I  rri  b  K 

3  0  3  0 

CRLL 

c 

di"' :R  =  f-,.  I=r.,Pl  X'  I  Ik  <<  r),R:nkx.  linkx 

3  0  4  0 

CRLL 

c 

rni.j  1  »  '  PO  ,  I  1  .  Pi  tik  ,  Is  nk  ,  Prnc  ,  I  rni;  ■ 

3  0  5  0 

Pnii:  ' 

I  > 

=  “k  ♦  P'rn i; 

3080 

Irr,.: 

I  ' 

=  -  k  »  I  m  c 

3  0  7  0 

R  m  d ' 

I  ' 

=  R  n  i  a '  I  '* 

3080 

I  m  d ' 

I  ' 

=  Imax  I  ■ 

3090 

NEXT 

I 

3  1  00 

R  m  a ' 

1 1 

'  =  R  rn  a '  7  ' 

3110 

I  ni  a ' 

1 1 

'  =  I  II  I  a '  7  ' 

3  1  20 

R  m  b  '■ 

1 1 

'  =  P  tti  b  '  7  ■ 

31  30 

I  m  b ' 

1 1 

'  =  I  rri  b  1  7  ' 

3140 

Rtrn:  ' 

1 1 

'  =  R  rri  ■:  1.  7  ' 

3150 

I  tri  i;  ' 

1 1 

’  =  I  m  '  7  ’ 

3  1  8  0 

R  m  d ' 

1 1 

'  =  R  m  d  '  7  , 

3  1  70 

I  md ' 

1 1 

>  =  I  rn  d '  7  ' 

3180 

R  rn  a ' 

10 

'  =  1 

3190 

I  rn  a  ■ 

1  0 

'=0 

3200 

Rmb  > 

1  0 

'  =  0 

3210 

I  rn  b ' 

10 

'=E.: 

3220 

Prnc  ’ 

1  0 

'=0 

32  30 

Im.:  ■ 

10 

'  =  0 

3240 

R  ft)  d ' 

10 

■  =  1 

3250 

I  rnd ' 

10 

=0 

B6 

::2d0 

R  n'l  -a  *.  9 

:;'  =  D 

Ji?0 

I  m  -a  ■  :  9 

■=0 

3230 

R  m  b  '■  9 

;■  =  0 

3230 

I  m  b  <  9 

i  =  0 

3300 

R  ni  c  <  9 

'=0 

3310 

I  me  '■  9 

=  0 

3  320 

R  m  d  9 

:■=£; 

3330 

I  md 9 

:>  =0 

5  3  4  0 

R  m  a  8 

'■  =Rmak  8  > 

3350 

I  mat  s 

■'  =  I  m  a  '■  8  ) 

33150 

R  ni  b ' .  S 

' =Rmb '  8 ‘ 

5  3  r  0 

I  Hi b ' .  S 

:■  =  I  rn  b  k  8  ' 

3  380 

Rmc  '  8 

>  =  R  ft  'i  *:  <  b  . ' 

3  3‘?0 

I  n'l  c  'v  o 

>  *  I  III  c  V  b  ‘ 

34  00 

R  n  i  d  '•  3 

■'  =  R  rn  d  k  8  . ' 

34  10 

I  m  d  '  8 

>  =  I  rn  d  ■  8  > 

34  20 

CALL 

Cdi  '.'kRjn'  5  .'  ,  lir-i'  5  ■  .  Rk  :  • 
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